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Abstract 

It was shown that at typical energies and luminosities of an e^e~ linear col- 
lider one can observe a significant suppression of fluctuations of the number of 
gluons in small angular regions. This property of jets is expected in e^e~ annihi- 
lation as well as in any collision producing high energy quarks which give rise to 
the QCD gluon cascade. The reduction of multiplicity fluctuations was obtained 
from perturbative QCD using bunching parameters and then cross checked with 
a Monte Carlo simulation. The parton-level Monte Carlo simulation confirms 
the prediction, while the hadron level illustrates that the effect survives the 
hadronization process. 
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1 Introduction 



A hadronic jet produced in high energy experiments can be viewed as a spray of 
strongly correlated particles created by a hard partonic scattering initiated at large 
momentum transfers A with A being a characteristic QCD scale) and subsequent 
parton cascade followed by a soft fragmentation process. At high energies, the dominant 
source of particle production inside jets is gluon splitting in the QCD cascade, where 
the presence of a gluon enhances the probability for emission of another gluon nearby 
in momentum space. This leads to inter-parton correlations and non-Poisson statistics 
for the multiplicity distributions in the restricted phase-space intervals where partons 
are counted. 

In the near future, experiments at the Large Hadron Collider (LHC) and, possibly, 
at an e^e~ linear collider (LC) will search for undiscovered Higgs bosons, new physics 
beyond the Standard Model as well as will focus on the precise measurements of known 
electroweak gauge bosons. Such studies will heavily rely on our understanding of the 
QCD multiple gluon radiation which is presently simulated with Monte Carlo (MC) 
models. These models take into account higher than first-order QCD effects in terms 
of the parton showers. This approach, however, has some shortcomings due to its 
approximate nature. Therefore, precise knowledge of gluon activity in jets together 
with reliable predictions of Monte Carlo models will require detailed and comprehensive 
studies of the jet structure. 

The internal structure of jets can be investigated using subjet multiplicities at some 
resolution scale Ucuti energy flows in specific kinematic variables with respect to the jet 
axis or calculating the jet shapes |T|. However, these studies do not address the issue 
of fluctuations of the number of separate hadrons in jets on event-to-event basis and 
do not resolve directly many-particle inclusive densities of the hadronic final state. 
To measure such fluctuations, one can calculate the probability distribution Pn{5) of 
observing n particles in a phase-space interval of size 5 inside a jet. Note that, in this 
case, we are interested not only in the average number {n{6)) of particles inside jets, 
but rather in the evolution of the shape of Pn{S) as a function of the size 6. This 
can be studied using the methods to be discussed below. Such studies have been done 
recently by the LEP and HERA experiments using various tools (see 0] for more 
references) . 

The semi-soft QCD physics, which determines multiplicities of the final-state hadrons 
inside jets, will become one of the important topics for investigation at LHC and LC. 
With an increase of the center-of-mass (cm.) energy of collisions, the multiplicities 
of produced hadrons inside jets will rise significantly. For example, for a LC with 
^/s = 500 GeV, the average number of hadrons in jets will increase by a factor two 
with respect to the LEPl experiments. This would allow detailed studies of the jet 
structure by counting separate particles inside jets and establishing relations between 
them in terms of correlations/fluctuations. Such a rise of the average hadron multiplic- 
ity, together with an increase in luminosity of the future experiments, will also allow 
the determination of the probability distribution Pn{S) for multiplicities higher than 
at LEP, assuming that the size 6 will be close to that of the previous experiments. 



1 



2 Methods 



One of the important problems in studying of multiparticle systems is to develop meth- 
ods which are sufficiently sensitive to the fluctuations of hadrons inside jets and make 
it possible an adequate comparison of theoretical predictions with experimental data. 
When fluctuations of separate particles produced inside jets are measured, it is conve- 
nient to transform the multiplicity distribution Pn{S) to the following observables: 



oo 



NFM: F,(5) = K5))^''5^-— -yP„(5), (1) 

n=q ^ '' 

CFM: K,{S) = F,{b) -Y, i/^'^^' (2) 

^-^ m\\q — m — 1 ! 

m=l ^ ' 

where the abbreviations denote the normalized factorial moments (NFM) P|, cumu- 
lant factorial moments (CFM) |]^ and bunching parameters (BP) [||. Note that the 
experimental deflnitions of the NFM and BP given in do not involve direct mea- 
surements of the Pn('^) and thus they differ from the theoretical formulae above. These 
three quantities measure deviations of the multiplicity distribution P„((5) from a Pois- 
son distribution P^°'*'^°"((5), since for this distribution Fq{S) = rjq{5) = 1, Kq{6) = 0. 
Note that such deviations are measured differently by these three methods |0. 

Uncorrelated particle production inside 6 leads to the Poisson statistics, thus devi- 
ations of the NFM, CFM and BP from the Poisson values mean correlations between 
particles and dynamical fluctuations. Let us note that, in a general case, uncorrelated 
phase-space production can have non-Poisson event-to-event fluctuations. In this case 
the NFM, CFM and BP are 5-independent constants which can differ from unity. 

In this paper we study the Pn{S) at asymptotically small 6, 

{n{6)) < 1, P„(5) > P„+i(5). (4) 

In this limit, the CFM are close to NFM and neither can resolve details of the high 
multiplicity tail of Pn(5), although both CFM and NFM are known to be sensitive 
to fluctuations of large multiplicities. Indeed, Eq. (^) means {n{6)) ~ Pi (5) and the 
deflnition of NFM can be approximated by 

i.e. the NFM are affected by the single-particle probability Pi (5). Thus even if the 
numbers of particles inside 6 follow a Poisson distribution for large multiplicities, Fq{6) 
are not unity if the Pi (5) is not determined by a Poisson value due to some kinematic 
reasons unrelated to multiple gluon emissions. 

This situation is avoided in the BP which measure the high multiplicities locally, 
near n — 1 multiplicity. (In fact, the BP are determined by the second derivative 
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Figure 1: A schematic representa- 
tion of the measurements of the 
probabihty distribution Pn{S,Q) 
inside jet. 



from lnP„((5) at sufficiently high n.) Therefore, even if Pn{S) for n < n' has a non- 
Poisson behavior, this does not affect the property rig{6) = 1 for events with particle 
multiplicities n > n' which obey the Poisson law. This important property of the BP 
will be used in the next section. 

3 Analytic QCD predictions for the BP 

There exist a number of the analytic QCD predictions for CFM and NFM of the 
parton multiplicity distributions obtained in the Double Leading Log Approximation 
(DLL A) [pHTll]. Here we discuss in more details the predictions derived in 0, since 
the others differ only by different treatment of non-leading perturbative contributions. 

The angular window considered here is defined for rings with the opening angle B 
around the outgoing-quark direction (see Fig. |I]). The phase-space interval where par- 
tons are counted is determined as half width 6 of the ring. Assuming sufficiently large 
energies, the probability Pn{S, 0) to find n partons inside the interval 26 is dominated 
by gluon bremsstrahlung off the outgoing quark. 

The analytic QCD predictions can be formulated in terms of the variable z, which 
is related to 6 as 

z = a-^ He/5), a = \n{Ee/A), (6) 

where E is the energy of the outgoing quark radiating soft gluons and A is the adjustable 
effective QCD scale. A decrease of the angular window 6 corresponds to an increase 
of the variable z. The maximum possible phase-space region {6 = 0) corresponds to 
z = 0. The variable z takes the maximum value z = 1 for 6 = A/E. 
The CFM derived in the DLLA satisfy the relation 

{n{6)yKg{z) = Ag exp [a{2 70(g) - z)] , (7) 

w, ~ VT^ (^q - (8) 

with 7o(Q) = >/2XT\q^{Q)Jt^ being the anomalous QCD dimension calculated at a 
scale Q ~ EQ, is the strong coupling constant and Ca = 3 is the gluon color factor. 
The functions Ag are determined by non-leading QCD contributions and reflect the 
values of CFM in the full phase space {z = 0). 

It is however clear from the previous section that the behavior of CFM at high 
orders q and large z (small S) is not strictly related to the structure of fluctuations 
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q\Pq{z) and thus the 



for multiphcities n > q. Similar to the NFM, their values are affected by the average 
number {n{z)) of counted particles which is dominated by Pi{z) for sufficiently large 
z. Therefore, to learn more about fluctuations without the bias from small-multiplicity 
events, one should calculate the BP. At large z, {n{z)YKq{z 
BP are 

A A -2 

and, assuming (H), one finally obtains 



)], 



(9) 



AqAg_2 



exp 



2a7o(Q)v/r^ln(l-z) 
g(g-l)(g-2) 



(10) 



where q > 2. It is remarkable that the expression in the exponent behaves as 1/g^ for 
large q. Therefore, the exponential factor for large q is very close to unity, i.e. 



AqAq_2 



q » 1. 



(11) 



This expression shows that the BP (0) for high orders q are close to those measured 
in the full phase space at z = 0, since they are only constructed from the functions 
Aq determined by non-leading QCD contributions. Generally, one expects the values 
of BP to be close to unity for 2; = 0, if the parton multiplicities in the opening angle 
G follow a distribution slightly broader than a Poisson distribution (i.e. a negative 
binomial distribution). 

A similar behavior of the BP is expected from alternative DLLA calculations. For 
instance, the NFM in the angular windows as defined above are given by fill 



Fq{z) = Fq{0) exp[z a (1 - Dq){q - 1)], 



Z},^2 7o(Q)- 



q + 1 fl-VT 



From (^ and the definition of the BP, one obtains 



(12) 



Vq{z) - ^?g(0) exp 



4a7o(g)(l- y/r^) 
g(g-l)(g-2) 

Fq{0)Fq.2{0) 



(13) 



(14) 



where again the exponential factor in ( |I3| ) is small at large q. Therefore, at large g, 
the BP are close to unity even if the angular phase-space region is very small. 

The behavior ( [TT| ) is opposite to that existing at small z (large 6). It is expected 
that the BP rise with increase of z, following a power-like behavior. Such a trend has 
been found recently by the L3 Collaboration for the rapidity distribution averaged 
over the all available phase space. The power-like behavior of the BP is also expected 
in simple cascade models [[T^ . 
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4 Monte Carlo tests 



The direct comparison of the analytic expectation discussed above with data is com- 
phcated by the fact that the calculation is performed in the DLLA which is know 
to contain a number of limitations. The most significant approximation is neglect of 
energy- momentum conservation in the gluon splittings (energy recoil effects). There- 
fore, it is important to compare the calculations with the MC parton shower which 
includes explicitly the recoil effect. In addition, comparisons with the hadron level of 
MC allow us to investigate the contributions from hadronization. 

To do this, we used the ARIADNE Monte Carlo program [|l^ for e~^e~ annihilation 
into hadrons. The hadronic events were generated at ^/s = 500 GeV with an integrated 
luminosity of 300 fb~^. The thrust axis was used to determine the jet direction. We 
used 9 = 25° with respect to the jet axis and the effective A = 0.15 to define the 
variable z. 

In addition to the hadronic sample, we generated a parton sample of MC events 
after switching off the hadronization phase. The perturbative phase was terminated by 
a cut-off in the transverse momentum of partons, Qq, which is chosen to be 300 MeV, 
to obtain the same average multiplicity of partons inside the opening angle as for 



hadrons. A similar method to analyze the parton correlations has been used in |TJ . 

Figure ^ shows the behavior of NFM for different orders g as a function of z, 
separately for hadrons and partons. Such a measurement, for q = 9, requires counting 
at least nine particles in small angular intervals. The smallest size of 6 for the present 
study is one degree {z ~ 0.5), which is close to that used at LEP and HERA 0. 

The factorial moments show a saturation at high z. The parton-level predictions 
systematically overestimate the moments measured using final-state hadrons. 

In addition to the MC predictions, the NFM obtained from the analytic DLLA 
calculations 0,|n| are also shown in Fig. 0. The major difference between these two 
analytic results is a different treatment of non-leading QCD contributions to the parton 
correlations. At present, the analytic calculations cannot give a satisfactory descrip- 
tion of the fluctuations, although they reproduce the rise of the moments. A similar 
conclusion has been made analysing the experimental data at LEP and HERA For 
the present comparisons with the MC model, we did not use the normalization constant 
Fq{0) (i.e. the NFM calculated at 2; = or B = Gq), in order to remove the theoretical 
uncertainty in the absolute values of Fg{z). Therefore, the theoretical curves shown in 
Fig. ^ have an arbitrary normalization, which, however, is not relevant for comparisons 
of the predicted z-dependence of the moments. 

Figure ^ shows the BP as a function of z for ARIADNE, separately for hadrons 
and for partons. The analytic estimates (|1^) and (|T3|) of the BP for q > 2 are also 
shown. The theoretical results can only be considered to be a realistic for large z values, 
because of the approximation Pn{z) ^ Pn+i{z) used to derive the BP from the CFM 



or NFM. It is important to note also that the analytic results (10) and (O) contain 



an unknown[] factor ?7g(0) ([T^), which is expected to be very close to unity. Therefore, 

^ It is possible to consider ln{riq(z)/riq{0)), in order to remove the theoretical uncertainty related 
to the full phase-space regions, similar to ln{Fq{z) / Fq(0)) used in For the BP, however, this leads 
to an unstable result due to different values of the BP at z = for partons and hadrons. Therefore, 
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calculating the analytic curves shown in Fig. a contribution of lnriq{0) is neglected. 
From this study, the following observations can be made: 

• High order BP for partons at large z show a clear tendency to approach the values 
close to those seen at z = 0. This is consistent with our analytic estimates in 
DLLA, suggesting a suppression of high-multiplicity fluctuations in small phase- 
space windows. 

• The hadron-level predictions show the same trend as partons. The values of BP 
for hadrons are closer to those for partons than for the NFM. This illustrates 
relative insensitivity of the BP to hadronization effects. The largest discrepancy 
between parton and hadron levels is seen at low z for r]2{z) which is affected by 
the probability Po{z) to find an event without particles. 

• For the analytic QCD calculations, there is a large difference between the two 
analytic DLLA QCD results, especially for the lowest order of BP. At large q 
and z, the analytic results and the MC predictions approach a common value 
(ln?7g(z) ~ 0). 

Thus we have established a very important observation which can shed light on the 
jet structure at small scales. According to the DLLA analytic calculations and our MC 
simulations, fluctuations of partons for sufficiently small angular windows are depleted. 

This study also illustrates that even for the relatively modest cm. energy and 
luminosity, the measurement of the NFM and BP can be possible up to g = 9, thus 
a LC will allow to observe the effect. So far the BP have been measured in e~^e~ 
annihilation up to g = 5, but the suppression of multi-gluon fluctuations in small 
rapidity regions has not been observed 0. 

5 High-multiplicity fluctuations in jets 

In this section we will summarize the observations made in this papers: 

The multiplicity fluctuations of particles for large phase-space regions {z ~ 0) are 
rather small, Fg{z = 0) ~ r]q{z = 0) ~ 1. Note that the fluctuations do not exhibit the 
Poisson structure since the BP and NFM show deviations of the NFM and BP from 
unity. 

With decrease of the angular window 6, the fluctuations start to emerge. At medium 
size of the window, z ~ 0.2 — 0.4, the fluctuations are maximal. 

For very small phase-space regions, the fluctuations are depleted and approach the 
values similar to those seen in the full phase space. This observation is only possible 
after the use of BP which are not affected by low-multiplicity part of the probability 
distributions, in contrast to the NFM and CFM. 

This peculiar feature of the fluctuations at low 6 can be understood in terms of 
strong restriction of the available phase space for the next generations of partons in 

such a normalization was not used in Figs. and pi 
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the cascade. One possible explanation of this saturation is to assume that particles 
counted inside 6 are produced by different parent partons and thus they are quasi- 
independent. An experimental verification of this phenomenon is rather important as 
it reveals the dynamics of the parton-shower system and can be directly related to the 
minimum transverse momenta Qo used in the perturbative calculations. 

The expected suppression of high-multiplicity fluctuations, however, seems to be 
different to that obtained in Refs. [^[I^) which also suggest a reduction of fluctuations 
in restricted pt intervals, or the existence of a critical angle 5^^^** beyond which the 
parton correlations are small and affected by the QCD cut-off. The latter effects should 
even be seen for the lowest order correlation functions and CFM, while the saturation 
discussed here is only expected in the tail of the multiplicity distribution measured 
with the BP. 

The fact that high-order BP are not sensitive to the hadronization dynamics for 
small 6 suggests applicability of the hypothesis of Local Parton-Hadron Duality ||T^ to 
fluctuations of high-multiplicity events. Therefore, the BP are the most suitable tool to 
study the high-multiplicity fluctuations and comparing experimental data with pertur- 
bative QCD predictions for partons. Some significant differences between fluctuations 
of partons and hadrons exist only for r)2, i.e. for low-multiplicity events. 
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Figure 2: The ARIADNE and analytic DLL A predictions for the NFM of different 
orders as a function of z defined in (|^) for e^e~ annihilation at = 500 GeV^. 
The lines show two different DLLA QCD predictions: (1) DLLA (solid lines); (2) 
DLL A (dashed lines). 
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Figure 3: The Monte Carlo and analytic DLLA predictions for the BP of different 
orders as a function of z defined in (^) for e~^e~ annihilation at ^/s = 500 GeV^. The 
lines show the BP from the analytic considerations: (1) DLLA of Eq. (|T0|) (solid lines); 
(2) DLLA of Eq. (|TB|) (dashed lines). Due to the approximation used to derive the BP, 
the analytic curves cannot be considered as a reliable prediction for small z. 
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